Fracture Surfaces as Multiscaling Graphs 
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Fracture paths in quasi-two-dimensional (2D) media (e.g thin layers of materials, paper) are an- 
alyzed as self-affine graphs h{x) of height /i as a function of length x. We show that these are 
multiscaling, in the sense that n*'' order moments of the height fluctuations across any distance £ 
scale with a characteristic exponent that depends nonlinearly on the order of the moment. Having 
demonstrated this, one rules out a widely held conjecture that fracture in 2D belongs to the univer- 
sality class of directed polymers in random media. In fact, 2D fracture does not belong to any of 
the known kinetic roughening models. The presence of multiscaling offers a stringent test for any 
theoretical model; we show that a recently introduced model of quasi-static fracture passes this test. 



Background: The pioneering experimental work de- 
scribed in Ref. P| drew attention to the fact that fracture 
surfaces are graphs in 2-1-1 (1+1) dimensions when the 
broken sample is three dimensional (two dimensional). 
This initial insight was followed up by a considerable 
number of works that focused on the scaling proper- 
ties of such graphs under affine transformations. In 1-1-1 
dimensions one denotes the graph as h{x) and considers 
the structure function S2{£), 



S2{e) = {\h{x + £) - h{x)\^) 



(1) 



where angular brackets denote an average over all x. The 
statement is then that this function is a homogeneous 
function of its arguments. 



(2) 



Close attention was paid to the numerical value of the 
scaling exponent C2/2 as observed in rupture lines in 
quasi two-dimensional materials (and in some numeri- 
cal simulations) H, Q, H, 1^ with the hope of assigning 
to fracture a "universality class" of one of the well stud- 
ied models of 1-1-1 dimensional kinetic roughening models 
Indeed, since in many such 1+1 measurements the 
exponent C2/2 was numerically close to 0.67, and also 
since it was proposed that fracture can be considered as 
a global minimization problem, some authors accepted a 
view that rupture lines map onto the model of directed 
polymers in random media i ^ ■ The latter model 
is consistent with an exactly soluble exponent C2/2 — 2/3 
0. In this Letter we show that this is not the case; more- 
over we advance reasons to believe that rupture lines do 
not map on any of the known 1+1 kinetic roughening 
models, but call for fresh thinking with new models in 
mind. One such new model will be shown to be a good 
candidate for this type of physics. 

Approach: We first argue in this Letter that the scal- 
ing properties of fracture lines are significantly richer 
than what can be read from Eq. |(2J). To this aim we 
consider the whole distribution function for height fluc- 
tuations, p{h{x + tj — h{x)j , and the associated higher 



order structure function Sn{£), 

Sn{e) = {\h(x + e) - h{x)n 



(3) 



For rupture lines the higher order structure function 
come each with its own exponent, 



5„(A^) ^ A«"5„(£) 



(4) 



where the scaling exponents ^„ are not simply related to 
the exponent (2 (i-e. Cn ^C2/2). In contrast, in di- 
rected polymers in random media the higher order struc- 
ture functions bring in no new information, and there 
Cn^nC2/2 M- 

Experimental example: An example that provides 
us with information of sufficient accuracy to establish 
the multiscaling characteristics is rupture lines in paper. 
The data acquired by Santucci et al. |lOj | was obtained 
in experiments where centrally notched sheets of fax pa- 
per were fractured by standard tensile testing machine. 
Four resulting crack profiles h(x) were digitized. Each 
digitization contained a few thousand points, where care 
was taken to insure that the smallest separation between 
points in x is larger than the typical fiber width; this 
is important to avoid the artificial introduction of over- 
hangs that destroy the graph property. 

Denoting Ah{£) = h{x + £) — h{x) we analyzed the data 
by boxing £ in logarithmic boxes, accumulating the data 
between 10° and 10°'^^ (the smallest box) and between 
10^-^^ and 10^-^ (the largest box). The m*'* box was 
considered as representing data for £ — io™^0-25. On the 
basis of this boxing we constructed the probability distri- 
bution function (pdf) P{Ah{£)) by combining data from 
all the four samples. Samples that exhibit marked trends 
(probably due to the finite size of the sample), were de- 
trended by subtracting the mean from each distribution. 
The computed pdf 's were then used to compute the mo- 
ments (O, and these in turn, once presented as log- log 
plots, yield the scaling exponents Cn- Such a typical log- 
log plot is shown in Fig. ^ exhibiting a typical scaling 
range of about 1.5 orders of magnitude. The resulting 
values of the scaling exponents are shown in Fig. [3 
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FIG. 1: A log-log plot of S3{£) as a function of £. The linear 
fit corresponds to a typical scaling range of about 1.5 orders 
of magnitude, with a slope of ^3 — 1.65. 



C =nH-n^X 



H= 0.64 ±0.02 
X = 0.026 ±0.002 




FIG. 2: The spectrum as a function of the moment order n 
for rupture lines in paper. The function is fitted to the form 
— nH — and the parameters H and A are given. The 
errors in the estimation of these parameters reflect both the 
variance between different samples and the fit quality. The 
linear plot nC,i is added to stress the non- linear nature of 



As a function of n these numbers can be fitted to the 
quadratic function C„ = nH — n?X with H = 0.64 ± 0.02 
and A = 0.026 ± 0.002 (a linear plot n^i is added for 
reference). The error bars quoted here reflect both the 
variance between different samples and the fit quality. 
The exponents were computed for n < 8 since for higher 
moments the discrete version of the integral 



|A/i(£)rP(A/i(£))dA/i(^) 



(5) 



did not converge. On the other hand, the convergence of 
the S*'* order moment is demonstrated in Fig. O 
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FIG. 3: An example of the convergence of the integral in Eq. 
for ^ = 10^-^^ and n = 8. 



The point to stress is that the scaling exponents C„ 
depend non-linearly on n; for the range of n values for 
which the moments converge, the exponents can be fit- 
ted to a quadratic function. It is well known from other 
areas of nonlinear physics, and turbulence in particular 
[Tl- 12. 13], that such phenomena of multiscaling are as- 
sociated with pdf's on different scales £ that cannot be 
collapsed by simple rescaling. In other words, in the ab- 
sence of multiscaling, there exists a single scaling expo- 
nent H with the help of which one can rescale the pdf's 
according to 



(6) 



where /(■) is a scaling function. In our case such rescal- 
ing does not result in data collapse. In Fig. 0]the nat- 
ural logarithm of P{Ah{i))£^ is plotted as a function 
of Ah{£)/e" for H = 0.64. Indeed, the data does not 
collapse onto a single curve. The fat tails of the proba- 
bility distribution functions at smaller scales is typical to 
multiscaling situations. 

The known kinetic roughening models, and in partic- 
ular directed polymers in random media, do not exhibit 
the multiscaling spectrum we have found, and thus we 
can infer that the scaling properties of rupture lines in 
1-1-1 dimension do not fall in the universality class of any 
one of the former. We propose that the long range elas- 
tic interaction that is typical to fracture is the origin of 
the positive correlations resulting in a non-trivial expo- 
nent (note that the quasi static fracture problem involves 
the solution of the bi-Laplacian equation) . To appreciate 
this important difference from standard kinetic roughen- 
ing models in 1-1-1 dimensions we turn now to a similar 
analysis of a theoretical model. 

Theoretical example: A recently published work 
(14l llSi] described a crack growth model aiming at un- 



3 



^ -4- 
< -5- 

-7- 
-8 
-9- 
-10 




-'15 -10 -5 5 

Ah{£)/£" 



10 



15 



FIG. 4: The natural logarithm of P{Ah{e))e" as a function 
of Ah{£)/£" for H = 0.64. The legend gives the scale £ for 
each plot. 



derstanding quasi-static fracture via voids formation and 
coalescence. The model consisted of two components: (i) 
the exact solution of the elasticity problem in the pres- 
ence of an arbitrarily shaped crack (ii) a growth law in 
which the evolution of the crack is controlled by the nu- 
cleation of voids at a finite distance R ahead of the crack 
tip. The first component involved the adaptation of the 
method of iterated conformal maps to the field of elastic- 
ity in complex geometrical domains |l6j . The second 
component incorporates an important physical insight 
about the role of a typical length scale R in the prob- 
lem. This length scale determines where voids nucleate, 
making the crack growth a succession of rapid growth 
steps interrupted by slower void nucleation events. This 
is a crucial factor that enables the crack to develop pos- 
itive correlations (i.e. a roughness exponent such that 
H > 1/2) [13. The succession of growth steps intro- 
duce variations in the crack geometry on the scale of 
the step size R. At a distance R ahead of the crack tip 
these geometric irregularities change the solution of the 
stress field and mediate positive correlations. When the 
growth takes place right at the crack tip Ellli, the 
crack appears locally straight; the geometric irregulari- 
ties are effectively screened and no positive correlations 
appear, i.e. one obtains H — 1/2. The origin of the finite 
length scale that is involved in the crack growth process 
can be the near tip non-linear physics (for example plas- 
ticity, see 0,01) or the scale of the quenched disorder 
in the system. The existence of a finite scale before the 
tip appears crucial for the development of a non-trivial 
roughness exponent. 

In [11 [13 it was shown that this particular growth 
model generates rupture lines with C2/2 « 0.66. It is 



thus interesting to examine whether this model exhibits 
the same kind of multiscaling that was found in the ex- 
perimental example above. Unfortunately, due to the 
significant computational cost of the iterated conformal 
maps technique the numerical investigation of the growth 
model had a limited number of realizations of a few hun- 
dreds growth steps. Due to the relative paucity of data 
the structure functions defined in Eq. Q would not con- 
verge well enough to provide reliable exponents. Instead, 
we use the max-min method which was proved reliable 
for the range of exponents under study and for the typical 
length of our cracks :20]. Thus, we define 

(7) 

The scaling exponents C„ are defined in analogous way 
to Eq. ^ by 



(8) 




FIG. 5: The spectrum („ as a function of the moment order 
n for rupture lines in the model of [T3. Il5|. The function is 
fitted to the form — nH — n^X and the parameters H and 
A are given. The errors in the estimation of these parameters 
reflect both the variance between different realizations and fit 
quality. The linear plot n(i is added for comparison. 

The resulting exponents Cn for the cracks generated 
by this model are shown in Fig. |S1 For the low or- 
ders moments (here we are limited by the paucity of 
data to n < 5) one again fits a quadratic function, with 
H ^0M± 0.02 and A = 0.023 ± 0.002. The errors in the 
estimation of these parameters reflect both the variance 
between different realizations and the fit quality. The n 
dependence of the exponents Cn and the values of the fit- 
ting parameters are in agreement with the experimental 
ones. Since there is nothing in the model that is specific 
for the physics of paper, it appears that multiscaling is a 
generic property of the fracture process, at least in l-fl 
dimensions. 

Summary: We examined an experimental example of 
rupture lines and a theoretical model of rupture in H-l 
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dimensions, and showed that both exhibit graphs whose 
scaUng properties appear similar, falling in a different 
class compared to standard models like directed poly- 
mers in random media. The scaling exponents of higher 
order structure functions associated with graphs of rup- 
ture lines depend non-linearly on the order of the mo- 
ments, while in directed polymers this is not the case. In 
fact, it appears that none of the standard known growth 
models in 1+1 dimensions is in the same class as the rup- 
ture lines discussed here. The analysis of the theoretical 
model appears encouraging, in the sense that it captures 
the multiscaling characteristic of the experimental exam- 
ple. We reiterate our proposition that the crucial aspect 
of the theoretical model that is responsible for the non- 
trivial scaling behavior is the existence of a finite scale R 
that does not exist in pure elasticity theory, or in models 
that treat the crack dynamics as a continuous process. 

Further experimental work is necessary to provide fur- 
ther examples of rupture lines in 1-1-1 dimensions. With 
more examples one would be able to determine the degree 
of universality of the scaling exponents, the dependence 
of the multiscaling exponents on issues like the isotropy 
of the cracking medium, the presence of plasticity and 
the protocol of the rupture experiments. 
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